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(D ' Abstract 
43 ■ 

We present a comparative study of the dynamical behaviour of topological sys- 
• i-H ! terns of recent interest, namely the non-Abelian Chern-Simons Higgs system and the 

Yang-Mills Chern-Simons Higgs system. By reducing the full field theories to tern- 
■ poral differential systems using the assumption of spatially homogeneous fields , we 

study the Lyapunov exponents for two types of initial conditions. We also examine 
in minute detail the behaviour of the Lyapunov spectra as a function of the various 
coupling parameters in the system. We compare and contrast our results with those 
for Abelian Higgs, Yang-Mills Higgs and Yang-Mills Chern-Simons systems which 
have been discussed by other authors recently. The role of the various terms in the 
Hamiltonians for such systems in determining the order-disorder transitions is em- 
phasized and shown to be counter-intuitive in the Yang-Mills Chern-Simons Higgs 
systems. 
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1 Introduction 



In recent times, the theory of dynamical systems has provided much insight into the 
origin of chaos in classical systems which were traditionally thought of as being completely 
deterministic. However, much of the progress has been mainly in the context of discrete 
mapping and those differential dynamical systems having low dimensional phase spaces. 

In the context of mathematical physics , there exist differential dynamical systems 
described by a large number of variables and therefore having phase spaces of rather 
large dimensions. Examples of such systems are the Yang-Mills system (YM), the Chern- 
Simons system (CS) and their various enlargements such as the Yang-Mills Higgs (YMH) 
, Yang-Mills Chern-Simons Higgs (YMCSH) and Chern-Simons Higgs (CSH) systems. 
These systems are treated as dynamical systems after being derived from the original 
highly nonlinear partial differential equations (PDE'S) through the assumption of spatial 
homogeneity which reduces the dependence of the dynamical variables on the three or 
four-dimensional space-time coordinates to a dependence purely on time. The systems 
then become temporal differential dynamical systems. 

Savvidy et.al jl| were the the first to demonstrate the chaotic nature of gaugetheories 
by treating them as differential dynamical systems. In fact, a detailed investigationhas 
been conducted in this context @ to classify the dynamical version of the pure YM field 
theory in terms of its ergodicity properties. 

While the condition for ergodic behaviour of a system is that the generic phase trajec- 
tory visits all regions of phase space given a sufficiently long time, and all phase-averages 
can be replaced by time-averages, the YM system exhibits stronger stochasticity prop- 
erties in phase space. Investigations revealthat it is definitely a mixing system, i.e., no 
time-averaging is required to achieve 'equilibrium'. In contrast to systems which are sim- 
ply ergodic its spectrum is continuous. Indeed, the studies seem to indicate that the YM 
dynamical system is a Kolmogorov K-system, which exhibits stronger mixing properties 
than those mentioned above. A connected neighbourhood of phase trajectories in this 
case, exhibits a positive average rate of exponential divergence or net positive Lyapunov 
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exponent. Equivalently, by a remarkable theorem a K system has positive Kolmogorov- 
Sinai (KS) entropy which is a measure of the degree of chaos, analogous to entropy as a 
measure of disorder in stastistical mechanics. In fact, while the KS entropy itself is not 
straightforward to measure, the fact that it is a sum of the positive Lyapunov exponents 
implies that it is an important concept in the classification of chaotic systems. 

Another aspect to the study of gauge theories as dynamical systems is the attempt 
to understand the ground (vacuum) state structure of quantum chromodynamics (QCD), 
and also its behaviour in extreme environments (such as high temperature). In this 
context, systems such as YMH, YMCS and YMCSH have been studied [Q. 

The addition of the CS term to various Abelian and non-Abelian gauge theories leads 
to novel features in general, as it is a topological term. Even in the complete PDE's, while 
a study of the YM case reveals interesting results in connection with the geometry and 
topology of four- dimensional manifolds, the related CS PDE's have yielded information 
about three-dimensional manifolds. The symplectic structure of CS theories differs in 
important ways from that of the Maxwell or YM gauge theories. In perturbative gauge 
theories the pure CS theories exhibit features that are absent in the Abelian gauge theory 
with both the Maxwell and CS terms. Delicate aspects relating to the infrared and 
ultraviolet behaviour and the regularization dependence in such perturbative theories 
||, the natural connection of quantized three-dimensional CS gauge theories with two- 
dimensional conformal field theories ||, and the effect of the YM term which acts as a 
singular perturbation when added to the SU(2) CS theory f7|, have all been extensively 
explored in the literature. 

Another aspect of interest in CS theories relates to its quantum mechanics. Non- 
perturbative quantum mechanical anomalies in these theories ||, infinite dimensional 
symmetry groups that arise in CS quantum mechanics 0], self-dual CS theories and 



extended supersymmetry ]T0| are a few areas in which distinct signatures of CS theories 
in sharp contrast to those of other gauge theories have been reported. 

In this paper we report on yet another aspect of CS theories and contrast them with 
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other gauge theories. This pertains to the chaotic nature of gauge theories mentioned 
earlier. For our purpose, the equations of motion are made to evolve only temporally, by 
suppressing the spatial dependence. While the resulting equations are like the continuum 
analogues of discrete maps (the latter being an area where extensive work has been done 
on their chaotic behaviour), they are also the dynamical version of the full gauge theory 
and hence represent one sector of the corresponding field theory. In the literature this is 
used as a convenient reduction to examine the integrability properties of the theory. This 
is because if this sector is proven to be non-integrable the corresponding field theory also 



will be chaotic [11]. 

In our recent papers, we have shown that the Abelian CSH system without a kinetic 
term is integrable, while the inclusion of a kinetic term, making it into a Maxwell Chern- 
Simons Higgs (MCSH) system (or Yang-Mills with a U(l) symmetry group) yielded a 



non-integrable system which admitted chaos [12|. The systems were also examined for 
the Painleve property. A numerical study of the Lyapunov exponents and phase space 
trajectories were carried out to show the existence of chaos in these systems. In ||13|| , 
the analysis was extended to the non-Abelian CSH and YMCSH system with an SU(2) 
symmetry group and both were found to be chaotic. 

In all these studies, one aspect that deserves more attention is the possibility of ob- 
serving the existence of a phase transition, i.e., is there a sharp order to chaos transition 
in the parameter space of these theories? Of course in these systems energy can also 



be used as a parameter since it is a conserved quantity. Indeed in ref.[|T4| , it has been 
claimed that an order-chaos transition is seen within a narrow range of energies. More 



recently, Kawabe JT5| has argued that the Abelian Higgs system (YMH with a U(l) sym- 
metry group) also shows an order-chaos transition for certain ratios of the two parameters 
in the theory. Our paper examines this aspect of chaos in YM systems by studying the 
non-Abelian CSH (NACSH) and the YMCSH systems with an SU(2) symmetry group. 
A comparative study is done to see the role of the kinetic term, the Higgs term and the 
CS term in the transition. 
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Some interesting features regarding the details of the phase transition from order to 
chaos in the dynamical analogues of both Abelian and non-Abelian gauge theories have 
been reported in the literature. In the context of Abelian Higgs theories Kawabe has 
reported transition from order to chaos within a certain range of the Higgs coupling 
constant and energy [15|. The onset of chaos is remarkably different qualitatively from 



the corresponding transition in the YMCS system where Giansanti and Simic [14| have 



reported the existence of an interesting 'fractal' structure in the phase transition region. 



Much earlier Savvidy et. al [16] have observed that the role of the Higgs is to order the 



Yang-Mills system and later extensive work on the YMH systems was conducted in ref 



17 ] . The picture that emerges therefore, is that the Higgs and the Chern-Simons term 
have distinct and different roles to play in the transition. It is therefore of importance to 
examine the effect of both terms on the Yang -Mills field in the YMCSH theory. As a 
primer to this, the competing effect of the 'oscillatory' behaviour of the CS term and the 
'stabilizing' role of the Higgs term in the CSH system is investigated in this paper. Later, 
we contrast this with the corresponding results in the YMCSH system. It is thus obvious 
that a K system like the Yang-Mills theory when coupled to CS and Higgs must show a 
rich and instructive behaviour in the understanding of regularity vs. chaos in Hamiltonian 
systems. A second aspect which emerges is related to an interesting question - will the 
Higgs stabilize any gauge invariant term involving only the vector potentials, independent 
of whether it is of the YM type or the CS type or is it necessary that the gauge field is 
only YM in nature? We attempt to answer these questions in this paper. 

2 The Dynamical Systems 

In this section we set up the two systems we shall be examining. Let us consider first the 
non-Abelian pure Chern-Simons Higgs system (i.e., without the kinetic term). From our 
earlier paper El, the Lagrangian for the non-Abelian (SU(2)) CSH (NACSH) system in 



2+1 dimensions in Minkowski space is given by: 



L = ^ w \f; v AI - y -f abc AlA b u A c a ] + D^ a D^ a - V(cf>) 



where 

F; u = d,A a v - d v Al + gfatcA^At, (2) 
f a bc are the structure constants of the SU(2) Lie algebra and 

D^ a = (d, - igT l A\)<\> a . (3) 

Here, T a are the generators of the SU(2) algebra, so that tr[T a Tb] = 5 a b- The equations of 
motion become: 



me^F^ = ig[D v ^T a (j> - tfT a D v <f>] (4) 

18V 

Yd4> 



D^{4>) = --— (5) 



Then, considering the real triplet representation for the Higgs field and the spatially 
homogeneous solutions diA a - = di<p = i,j=l,2 and the gauge choice Aq = we get for 
the v = component of eqn. [4] 

mA\ x A 2 = —0 x cf) (6) 

which is just the Gauss' law constraint. The remaining equations of motion for the vector 
field are: 

A\ = 9 -{A 2 ft - U2 ■ $) (7) 
m 

a\ = ——{M(f 2 - 4>Ai ■ $) (8) 

m 

The equation of motion for the Higgs field is: 

$= -g 2 [(Aj + 4)0 - (A, ■ $Ai + M ■ $M - (9) 

From the equations of motion for the vector fields it is easily seen that A\ + A\ is a 
constant of the motion. Throughout this paper, we work with the potential 

V(0) = ^(0 2 -„ 2 ) 2 . (10) 

The NACSH system described by the above equations of motion has three parameters. 
For comparison with the work of Kawabe | 15| , we shall scale the variables such that we 
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are left with only one parameter. The following scaling of variables, 

A 1 — ► gA 1 
A 2 — ► gA 2 



m 
9 



reduces the equations of motion to: 

M = [A 2 fi - ${A 2 ■ $)\ (11) 
A 2 = -[A^ (12) 
4> = -[(Al + A 2 2 )$-(A 1 -$A 1 + A 2 -$A 2 )]-^ 2 -v 2 ) (13) 

with k — ^ . In the rest of the paper, we shall set the scaled v to be one without 
loss of generality. The scaled Lagrangian is given by 

L = (X ■ A 2 - A 2 ■ Ax) + 2 - [{A\ + 4)0 2 - {M ■ 0) 2 - (A 2 • 0) 2 ] - ^(0 2 - l) 2 (14) 

The corresponding energy function is easily seen to be 



-?2 



E = + [{A\ + A|)0 2 - (A*l • 0) 2 - (A 2 • 0) 2 ] + ^(0 2 - l) 2 . (15) 

These are the equations governing the NACSH dynamical system. Since we wish to 
compare the results for the NACSH system with those of the YMCSH system we proceed 
to set up the YMCSH dynamical equations. The Lagrangian is given by 



L = -^F^ + ^iF^-p^AlAlA^} 

+D^ a Dy a - V(4>) (16) 



The equations of motion are 

D,F^ a + me^F^ = ig[D l '^T a <f> - <pT a D"4>\ (17) 

dV 1 dV 
<90* 2 d<p a 
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The v — component gives the Gauss' law constraint which in this case is 

-(Ai x Ai + A 2 x A 2 + 2mA 1 x A 2 ) = -20 x 0. (19) 

Once again, we choose the gauge A = and consider the spatially homogeneous case ; 
then the equations of motion for the gauge fields become: 

A 1 + 2m A 2 + 2g 2 (A 1 <f - $A X ■ 0) 

+g 2 (A 1 A 2 ■ A 2 - A 2 A X ■ A 2 ) = (20) 
A 2 - 2m A 1 + 2g 2 (A 2 (f 2 - $A 2 ■ 0) 
+g 2 (A 2 A 1 ■ A 1 - A 1 A l ■ A 2 ) = 0. (21) 

From these equations it is easy to see that 

A 2 ■ A 1 - A x ■ A 2 + m{A\ + A 2 2 ) 

is a constant of the motion. The equation of motion for the three component Higgs field 
becomes 

$= -g 2 [{A 2 + 4)0 - {A, ■ Ui + M ■ }A 2 )\ - (22) 

1 (70 

In this case we have three second order differential equations for each vector field. The 

Lagrangian which leads to these equations of motion is: 

\ '-> '-> ' ' 

L = -(A 2 + A 2 ) + mil, ■ A 2 - A 2 ■ A,) + 

-g 2 [\{A 2 A 2 - (A, ■ A 2 ) 2 ) + (A 2 + A 2 )fi - (A, ■ 0) 2 - (A 2 • 0) 2 ] - V{4>) (23) 

Using the same rescaling as for NACSH we have the equations of motion 

-A\ + 2A 2 + 2(A 1 0~' 2 - 4>A X ■ 0) 
m 

+—(A-yA 2 ■ A 2 - A 2 Ai ■ A 2 ) — (24) 
m 

—A 2 - 2A X + 2{A 2 $ 2 - IA 2 ■ 0) 
m 

+-(A 2 A\ ■ A ± - A^ ■ A 2 ) = (25) 
m 



and 



$ = -[(if + A 2 )$ - (A, ■ + M ■ 0A 2 ] - ^0(0 2 - 1). (26) 



It is interesting to note that while in the NACSH system the Yang-Mills parameter g, the 
Higgs parameter A and the Chern-Simons parameter m could all be combined into the 
parameter k, this is not possible for the YMCSH system where we are left with both k 
and m appearing explicitly. The scaled Lagrangian in this case is given by 

L = ^(if + A\) + (lx • A 2 - X 2 ■ A\) + / 

- (Ai • Mf] + [(if + 4)0 2 - {M ■ 0) 2 - (A • 0) 2 ] - - l) 2 (27) 

with an energy function 

E = i^(^+Al)+t +^[Aaia^ 

(28) 

This completes the description of the dynamical systems which we shall be studying. In 
the next section, we describe the numerical analysis that we have carried out. 

3 Numerical Analysis 

It is clear from the above dynamical equations that the phase space of these systems are 
large. One of the traditional ways for examining phase spaces to determine chaotic vs. 
regular behaviour has been to use the Poincare sections where one examines the points 
mapped out on a plane surface in phase space as the trajectory crosses it. While this is 
straightforward for dynamical systems whose phase space dimensionality does not exceed 
four, it is difficult to interpret it in the systems we are dealing with. 

Instead, we examine the variation of the maximal Lyapunov exponent as the two 
NACSH parameters energy and k are varied. This clearly shows us regions of regular 
behaviour (where the exponent goes to zero) and regions of chaotic behaviour (where 
the exponent is positive). These calculations were carried out for a wide range of initial 
conditions. The initial conditions that were chosen were in turn dictated by the dynamical 
systems themselves. Being derived from the equations of motion the field variables are 
required to satisfy the Gauss' law constraint. Since this constraint must be preserved 
during the time evolution of the dynamical system, it is sufficient to ensure their validity 
via the initial conditions. 
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For the NACSH system the following forms were chosen as the initial conditions: 




A, 







' -x \ 


1 


/ X 


X 


H 




d> = - 
r 2 


X 


{ o J 











(29) 



For YMCSH this is supplemented with 

Ai 
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v o y 




\0 



(30) 



x was then varied to obtain a range of initial conditions for suitable energies of interest. 
In figs. 1-4 we show the behaviour of the maximal Lyapunov exponent as a function of 
energy for k = 0, .5, 1, 5. For k = (i.e., absence of the Higgs potential) we see that the 
system is mostly chaotic with a window of regularity for 7 < E < 9. Increasing the value 
of k , we find a transitional region of regular to chaotic behaviour at small energies, in 
contrast to mostly oscillatory behaviour at higher energies. For k — 1, more transitions 
from order to chaos appear at larger energies. Much more dramatic behaviour is seen for 
large k (k — 5), where oscillatory behaviour manifests itself for E > 3. Thus we see that 
the effect of the topological term (large m) is to produce a regular oscillatory behaviour 
in the dynamics of the system. 

To illustrate the regular and chaotic behaviour on the trajectories of this system we 
exhibit in figs. 5 and 6, phase plots corresponding to the energies at which the Lyapunov 
exponent shows regular and chaotic behaviour. These correspond to k — 1 and energies 
10.28 and 20.607. 

Giansanti and Simic |14| report a fractal-like structure of order-chaos transitions in 
YMCS systems. In our particular case, when the Yang- Mills field is absent and the Higgs 
field is present we see no such fractal behaviour in the region of phase space that we have 
examined. Therefore, this suggests that the quartic coupling that arises from the inclusion 
of the kinetic YM term may be responsible for the observed fractal structure. This rich 
phase space structure of the CSH, YMCS and YMCSH systems clearly needs further 
exploration. The set of trajectories examined by Giansanti and Simic do not correspond 
to the ansatz that we have chosen and therefore we have explored different regions of 
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phase space. Thus our results are complementary to those obtained by Giansanti and 
Simic. 

We have examined the NACSH system for another ansatz which we may call the two-v 
variable ansatz: 

Al -{1) f -[l) m 

For YMCSH this is supplemented with 















R 






(32) 
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V o y 





This allows us to compare our results with that of Kawabe for the Abelian Higgs theory 
[T5|j , where motion is completely bounded for Q = < 1, while it is unbounded for 
Q > 1. In contrast, for the CSH system that we are studying , for all values of Q = 4E/ 
k and k = -Ay, the motion is bounded. Indeed we find that for large E and large k the 

_,2->2-> -> -> -> -> -> 

contour levels of the function W = [(A 1 + A 2 )4> 2 - (A 1 ■ (p) 2 - (A 2 ■ (p) 2 } + f (4> 2 - l) 2 show 
extremely restricted domains for the dynamics. Fig. 7 represents the potential contours 
for k = 10 for energies in the range {0.25, 2} while fig. 8 gives the contours for the energy 
range {2, 10}. There is a dramatic change in the available phase space for this particular 
two- variable ansatz. This clearly suggests that the non- Abelian nature of the CS term is 
already contributing to a significant change in the dynamical structure. 

Evidence for the transition from regularity to chaos is seen in fig. 9, where we show the 
fraction of the phase space that is regular. This is obtained by calculating the maximal 
Lyapunov exponent for k = 1 for various initial conditions for energies ranging from 1 to 
10. The cutoff on the exponent for regularity was taken to be 0.01. A simple count on 
the exponents falling below this value out of a hundred initial conditions for each energy 
was carried out. The figure clearly reveals that for small k , the NACSH system becomes 
chaotic more or less monotonically as the energy is increased. However for large k there 
is no such simple behaviour . A fuller discussion of this is given in the next section. 

This transition between regular and chaotic behaviour is also seen in the phase space 
trajectories shown in figs. 10 and 11. In fig. 10 we plot 02 vs. <pi for k — 1, energy=10 
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and an initial condition where the maximal Lyapunov exponent is almost zero. In fig. 11 
we show the phase plot for the same parameter values for an initial condition which gives 
a large maximal Lyapunov exponent. 

We now proceed to the numerical investigation of the YMCSH systems. Fig. 12 shows 
the variation of the maximal Lyapunov exponent as a function of energy for k = 1, 10 and 
energy=10, with m=l. The graphs clearly show that for large k (where either the Higgs 
coupling A is large or the YM coupling g is small) the system exhibits more regularity for 
low energies. 



This is to be contrasted with earlier work |T|],[ 15 1, where it was reported that in a YMH 
system as A increases the system becomes more regular regardless of the energy regime 
being investigated. Here, we see that for the YMCSH system such transitions to regularity 
are seen for small energies as k increases while the Lyapunov exponent increases almost 
linearly with energy for the large energy regime. While fig. 12 was obtained for the one- 
variable ansatz, fig. 13 shows the maximal Lyapunov exponent as a function of the initial 
variable x for the two-variable ansatz, for fixed values of energy and k. Once again we see 
regions of regularity for small energies and chaotic behaviour for large energies irrespective 
of the value of k. Evidence for regularity and chaos for the two-variable ansatz is given in 
figs. 13 and 14. In fig. 14 we see that for k = 5 and energy=l, the phase space trajectory is 
highly regular and quasi-periodic. Fig. 15 shows a region for the same parameter values, 
where the phase space is chaotic. 

All the calculations for the NACSH system were carried out using a straightforward 
Runge-Kutta fourth order routine with care being taken to preserve the constants of 
motion to an accuracy of one part in 10 5 . However, the calculations for the YMCSH 
system required an adaptive step-size Runge-Kutta routine to ensure energy conservation 
to the same degree of accuracy as in the NACSH system. 
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4 Results and Discussion 



From the numerical studies undertaken by us, certain very interesting features emerge, 
not only regarding the richness of the phase space corresponding to gauge theories, but 
also with respect to some 'counter-intuitive' phenomena that occur in these systems. 



One feature that seems to be common to both the YMH [16] and CSH systems, from 
an investigation of the potential contours corresponding to various initial conditions, is 
the boundedness of phase space. This is in contrast to the situation that prevails in the 
Abelian Higgs system [[HJ. The boundedness observed is better understood when we 
realize that the YMCS dynamical system can be used to describe particle motion in a 



magnetic field a physical situation which allows only for bounded motion. 

Apart from the crucial role played by the ansatz chosen by us, in this matter, the non- 
Abelian nature of the gauge term also plays an important part. In this context, we realize 
that the boundedness observed is quite independent of the details of the non-Abelian 
coupling, i.e., whether it is of the YM type or of the trilinear CS type. 

We recall that the YM system is a K-system || and that the CS term creates more 
regular windows when added to the YM system [14| . This accounts to a large measure, for 
the fact that for some k and m values, the maximal Lyapunov exponents in the YMCSH 
case are larger than the corresponding ones for the CSH system. 

However, a striking feature that emerges in our studies that is counter-intuitive is that, 
in general, it is not true that increase in k 'regularizes' the gauge term at all energies. 
Recall that our calculations are made in terms of the scaled parameter k. An increase in 
k could either be due to an increase in the Higgs coupling A or a decrease in the gauge 
coupling g, for a fixed m. Whereas in the former case the regularity is expected to increase 
, in the latter case, it is not established that for all small non-zero g more regular islands 
appear. In fact, as borne out by figs. 12,13 and 14, an increase in k produces more 
regularity only for small values of energy. This can be understood better if we 
realize that it is not just the value of k that determines the appearance of regular islands, 
but also the available phase space as well. This is clear from the energy contours shown 
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in figs. 7 and 8. With more accessible regions in phase space, it is but natural that more 
randomness sets in and the KAM tori structure gets broken. 

Another aspect is that as the energy increases, the maximal Lyapunov exponent in- 
creases in magnitude in the YMCSH system. This shows that the YM terms takes over 
for large energies and the CS term produces the 'oscillatory effect'. The effect of the CS 
term is reminiscent of the 'fractal' structure observed in YMCS systems, where in various 
energy windows, order-chaos-order transitions are observed 0]. 

An important feature that we see is that while in the pure CSH system large windows 
of regularity exist ( as seen for the range of energies shown in figs. 1-4), in YMCSH this 
does not happen. This perhaps can be traced back to the interplay between the oscillatory 
effect of the CS term, the regularizing effect of the Higgs term and the completely chaotic 
nature of the YM term. Further , with increase in k , transitions from regularity to chaos 
take place in the pure CSH case at relatively higher energies when compared with the low 
energy values corresponding to order-chaos transitions in YMCSH. 

The final picture which emerges bears out the fact that in a complex dynamical system 
with a large phase space (in contrast to the wide class of Hamiltonian systems with two 
degrees of freedom) curious interplay between different coupling constants and the rich 
structure of phase space itself can lead to novel results- some of them quite counter- 
intuitive and surprising. A detailed examination of dynamical systems which emerge 
from field theoretic systems are thus vital for an understanding of Hamiltonian systems 
with a large number of degrees of freedom. In turn it is also an important primer in the 
understanding of non-Abelian field theories themselves. 
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6 Figure Captions 

Fig.l Variation of the maximal Lyapunov exponent with energy for NACSH for k — 0. 
Fig. 2 Variation of the maximal Lyapunov exponent with energy for NACSH for k = 0.5. 
Fig. 3 Variation of the maximal Lyapunov exponent with energy for NACSH for k = 1.0. 
Fig. 4 Variation of the maximal Lyapunov exponent with energy for NACSH for k = 5.0. 
Fig. 5 Phase plot of favs.ifii showing regular behaviour for NACSH for the one-variable 
ansatz for k = 1, E=10.2088 and x=1.381. 

Fig. 6 Phase plot of fovs.fa showing chaotic behaviour for NACSH for the one-variable 
ansatz for k = 1, E=20.6071 and x=1.64. 

Fig. 7 Contour plot for the potential corresponding to NACSH for k — 10 and energy 
range {.25,2} 

Fig. 8 Contour plot for the potential corresponding to NACSH for k — 10 and energy 
range {2, 10} 

Fig. 9 Phase space fraction of regularity vs. energy for k — 1 for NACSH 

Fig. 10 Phase plot of 2 ^s.0i showing regular behaviour for NACSH for the two- variable 

ansatz for k = 1, E=10, x=.05 and y=2.598405 

Fig. 11 Phase plot of favs.fa showing chaotic behaviour for NACSH for the two-variable 
ansatz for k = 1, E=10, x=.9 and y=1.908953. 

Fig. 12 Comparison of maximal Lyapunov exponent vs. energy for YMCSH for the one- 
variable ansatz for k — 1, m=l (* represents the curve) with k = 10, m=l (+ represents 
the curve). 

Fig. 13 Comparison of maximal Lyapunov exponent vs. x for YMCSH for the two- variable 
ansatz for k — 5, E=l,m=l (o represents the curve) with k = 10,E=5, m=l (+ represents 
the curve) and k — 1, E=10, m=l (* represents the curve) 

Fig. 14 Phase plot of fovs.cpi showing regular behaviour for YMCSH for the two-variable 
ansatz for k = 5, E=l, x=.725 and y=0.742085. 

Fig. 15 Phase plot of 02^s-0i showing chaotic behaviour for YMCSH for the two-variable 
ansatz for k = 5, E=l, x=.9125 and y=0.571046. 



18 



19 



